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ABSTRACT 

We analyze our accurate kinematical data for the old clusters in the inner regions of M31. These 
velocities are based on high S/N Hectospec data (Caldwell et al 2010). The data are well suited 
for analysis of M31's inner regions because we took particular care to correct for contamination by 
unresolved field stars from the disk and bulge in the fibers. The metal poor clusters show kinematics 
which are compatible with a pressure-supported spheroid. The kinematics of metal-rich clusters, 
however, argue for a disk population. In particular the innermost region (inside 2 kpc) shows the 
kinematics of the X2 family of bar periodic orbits, arguing for the existence of an inner Lindblad 
resonance in M31. 

Subject headings: catalogs - galaxies: individual (M31) - galaxies: star clusters - globular clusters: 
general - star clusters: general 



1. INTRODUCTION 

Globular clusters can provide simultaneous estimates 
of velocity, metallicity and age: a powerful trio with 
which to study the history of a galaxy. They are par- 
ticularly helpful to complement integrated light studies, 
which average over all stellar populations along a line 
of sight. In this paper we discuss the kinematics of old 
clusters projected on the inner 10 kpc of M31. Roughly 
one third of our sample of over 300 old clusters in M31 
(presented in Caldwell et al. 2009 and 2010, Papers 1 & 
2 hereafter) are located within 3 kpc of its center. Be- 
cause of our careful treatment of the effect of field star 
contamination from the bright bulge and inner disk re- 
gion in our fibers, our dataset is particularly well suited 
for study of the central regions of M3L 

Early work on bulge kinematics (jDavies "etld1[l98l 
showed that bulges resemble low-luminosity ellipticals in 
being kincmatically hot with a high degree of rotational 
support (V/a ^ 1). Later studies, however, showed that 



bulges are more complex and that an important distinc- 
tion must be made between classical i? 1 / 4 bulges - which 
are kinematically hot and formed rapidly from mergers 
and collapses - and bulges formed via secular evolution 
of disks, which have a lower Sersic index (Kormendy 
k, Kennicutt 2004). In this second category Athanas- 
soula ( 2005) distinguished the boxy/peanut bulges - 
which are parts of bars seen edge-on - and the disk-like 
bulges, which have a disk shape. Boxy/peanut bulges 
can be distinguished in near- edge-on galaxies from pho- 
tometry or via kinematics (eg Kuiiken fc M cmficldl ll995l ; 



iBureau fc Athanassoula|[i999ir 

Evidence from isophotal twists and kinematics was 
used to argue that M31 might have a triaxial bulge 
or a bar (Lindblad 1956, Stark 1977, Stark & Binncy 
1994). More recently, Athanassoula & Beaton (2006) and 
Beaton et al. (2007), using deep 2MASS observations, 
considerably strengthened the case for a bar and sug- 
gested that M31 also has a centrally concentrated classi- 
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cal bulge, which dominates the light in the inner 200 pc. 
We note here that this is a considerably smaller and less 
dominant class ical bulge than t he on e suggested by pre- 
vious authors: Ide Vaucouleursl (11958ft found an effec tive 
radius of 3.5 kpc. and iWalterbos fc Kennicuttl (|1988| ) de- 
rived an effective radius of 2 kpc and found that the bulge 
contributed 40% of the light of the galaxy. Our globular 
cluster kinematical data allow us to further explore this 
shift in our view of M31's bulge, since we have [Fe/H] 
and velocity measurements for 98 old clusters projected 
within 3 kpc of M31's center. 

We assume a distance of 770 kpc throughout 
(|Freedman fc Madord |1990[ ) and a PA of 37.7 degrees. 
The XY coordinate system we use in this paper has units 
of kpc, with positive X along the major axis towards the 
NE. 

2. CLUSTER KINEMATICS 

Paper 2 presented [Fc/H], age and v elocity measure- 
ments based on high S/N Hectospec ([Fabricant et al.l 
120051 ) spectra (a median S/N of 75 per A) for over 300 
M31 clusters with ages greater than 6 Gyr. (In fact the 
great majority of these clusters have ages greater than 10 
Gyr.) Here we discuss the old clusters from this paper 
which arc within 2 kpc of M31's major axis. Repeat Hec- 
tospec observations showed a median velocity error of 6 
km s _1 . Our study contains 17 entirely new cluster ve- 
locities and is the first fiber study to use offset exposures 
near each cluster in the bright inner regio ns to correct for 
the co ntamination from field stars there. iCaldwell et al.l 
(|2010D showed that ignoring this effect can lead to ve- 
locity errors of more than 100 km s _1 . In the small 
number of cases where our velocities differed significantly 
from the Hcctochclle velocities of Strader and Caldwell 
(2010, in preparation), we have used the more accurate 
Hectochelle data. Our [Fc/H] values are in good a gree- 
ment with the recen t results of iBeaslev et al.l (|2005f ) and 
iColucci et all (|2009D and the HST color-magnitude de- 
rived values. We found that the old cluster mctallicity 
distribution was neither unimodal nor simply bimodal, 
showing a median [Fe/H] around -1.0 and possible peaks 
at [Fe/H] = -0.3, -0.8 and -1.4. 

Previous work on M31 globulars suggested a larger 
systemic rotation for the metal-rich er clusters (eg 
iHuchra et ai][l99lt iBaTmbv et al.l [2000h . Since all but 
one of the clusters with [Fe/H] > —0.4 are projected on 
the inner disk (less than 1.5 kpc from the major axis) 
we first explore connections between the metal-rich clus- 
ters and M31's disk. In the Milky Way, the metal-rich 
globular clu sters are likely associated with its inner disk: 
Zinnl (119851) connect ed these clusters with the thick disk, 



Minniti et all (|1996D with the bulge. Since more recent 



work has shown that the Milky Way bulge is domi nated 
by a bar (jWeiland et al.l 11994b iBinnev et al.lll997D . the 
most metal-rich globulars in the Milky Way are all likely 
to be connected to its disk in some way: either as a bar 
distortion or a thicker component. 

To compare our old clusters with M31's disk, we use 
the fibers which were designed to measure the contami- 
nation from M31's disk and bulge: each fiber's measured 
velocity will be the mean velocity of all the stars along 
that line of sight. Figure [1] compares cluster velocities 
(with different colors for clusters with different metallic- 
itics: red for most metal rich through blue for the most 



metal poor) with these mean velocity estimates, shown in 
black. The top panel shows clusters with [Fc/H] > —0.6; 
the bottom panel shows more metal poor clusters. The 
difference between their kinematic behavior is stunning. 
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Fig. 1. — Velocity of clusters projected less than 2 kpc from 
the major axis. The upper panel shows the metal-rich clusters 
([Fc/H]>-0.6) and the lower panel the more metal-poor clusters 
which dominate M31's old clusters. In both panels, black symbols 
show the mean velocity of stars at that position, integrated along 
the line of sight. Open symbols denote ages less than 10 Gyr, 
closed symbols greater than 10 Gyr (note that we are unable to 
measure ages for clusters with [Fe/H] less than —1.0, and we use 
closed symbols for these clusters) . The vertical grey lines show the 
end of the th i ck ba r (das hed lines) and the thin bar (so lid), from 
IBeaton et al.l 1120071 1 and lAthanassoula &~B caton 120061) . and the 
solid black line is the rotation curve from IKentl 119891 V Note that 
measurements of the dimensions of the thin and particularly the 
thick bar are approximate only. 



Metal poor clusters (lower panel) show little sign of 
rotation and occupy the four quadrants of the plot simi- 
larly. On the other hand, the metal-rich clusters (upper 
panel) show a distinct and quite cold kinematical signa- 
ture. There arc almost no clusters in the forbidden quad- 
rants (occupancy here corresponds to rotation in the op- 
posite direction to the disk) and most of those more than 
2 kpc from the center (ie \X\ > 2) have velocities which 
closely follow the disk velocity at that position. However, 
in the inner 2 kpc the signature differs from the usual 
one for a disk composed of stars on near-circular orbits. 
Although all except one cluster occupy the same quad- 
rant as the disk, thus respecting the same direction of 
rotation, their velocities can deviate from the local mean 
velocity of the integrated light by up to 350 km s _1 (re- 
call that M31's rotation velocity is 250 km s _1 ). We note 
that very high velocities are also obse rved in the HI gas 
in this region (Brinks fc Shanelll984D . In the following 
section we will describe expectations for the kinematics 
of thin disk, bar and classical bulge objects, and show 
that this signature is expected for bar orbits. 
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2.1. Kinematics: expectation from disk, bar and bulge 

Thin disks in galaxies have "cold" kinematics domi- 
nated by rotation and show a low velocity dispersion. 
We showed in Paper 1 (see Figure 13) that the young 
M31 clusters (with ages less than 2 Gyr) have such kine- 
matics: the young clusters all follow the same narrow 
locus in position vs velocity. We also showed (see Figure 
12) the mean velocity field across the face of the disk, ob- 
tained from our "sky" fibers. The mean velocity changes 
smoothly and slowly as we look from the receding side of 
the disk through the center to the approaching side, as 
expected for a thin disk. 

It is particularly simple to follow the kinematic signa- 
ture of a cold, thin disk by examining velocities of objects 
seen close to the major axis. In a galaxy close to edge-on 
such as M31, such a star in a circular orbit will have all its 
velocity in the line of sight, giving a clean measure of V^, 
the azimuthal component, from the line-of-sight velocity. 
For disk stars observed at larger distances from the ma- 
jor axis, less of their azimuthal velocity will be projected 
onto the line of sight and so the change in mean velocity 
from one side of the disk to the other will be smaller. 

Most orbits in bars follow the two main famil ies of 
closed periodic orbits ([Binnev and Tremainel [2008T) : the 
X\ orbits, which are aligned along the long axis of the bar 
(close to the major axis in M31, see Beaton et al (2007)) 
and xi orbits which are aligned along its short axis (close 
to M31's minor axis). For X2 orbits, velocities can reach 
very high values close to the center of the galaxy. This 
is due to the fact that they are observed near-end-on, so 
that the line-of-si ght component is nearly along the orbit 
at its pe ricenter ([Bureau fc: Athanassoula l ll999D . Figure 
[5] (from IBinnev et al.1 119911) illustrates the spatial and 
velocity signatures of X\ and x^ orbits. It can be seen 
that in this example, the x^ orbits reach velocities much 
higher than the circular velocity. (A similar position- 
velocity diagram for a M31-like system can be seen in the 
midd le top panel of Figure 11 of lAthanassoula fc Beaton] 
[2003 ). 
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Fig. 2. — Illustration of regions occupied both spatially and 
kinematically by x \ and X2 orbits in a barred potential, from 
IBinnev et al.l 1 119911) . Left panel shows their spatial location in a 
face-on view, while the right panel shows longitude-velocity plots. 
x\ orbits are aligned along the bar major axis and shown with 
solid lines, while X2 orbits align perpendicular to the bar major 
axis and are shown with dotted lines. Note that in this exam- 
ple, the X2 orbits can reach velocities significantly higher than the 
circular velocity, which is v=l in this model. 

Lastly, we would expect any classical bulge component 



to show V I a ~ 1: some rotational sup port but a r oughly 
equivalent amount of random motion. iKentl ([1989) fitted 
the M31 bulge using an oblate rotator model with major- 
axis velocity of around 90 km s" 1 and velocity dispersion 
of 130 km s _1 at 1.5 kpc from the center. 

3. DISCUSSION 

We saw in Figure Q] that the kinematics of old M31 
clusters with [Fc/H] > —0.6 in its innermost region show 
the distinctive behavior of objects on X2 orbits in M31's 
bar. This is in very good agreement with orbital struc- 
ture in bars since the xi orbits are always confined to 
the innermost regions, in the region interior to the inner 
Lindblad resonance (ILR) . Most of the rest of the metal- 
rich clusters have orbits consistent with disk objects. 

We see little or no indication in the kinematics in the 
upper panel of Figure [T] for a kine matica l ly ho t popu- 
lation such as the classical bulge of [Kent (1989). How- 
ever, we note again that the classical bulge identified by 
iBeaton et al.l ([20071 ) was quite small, only dominating 
the inner 200 pc. We have only one cluster within 200 
pc of M31's center in our sample, so we cannot probe the 
kinematics of this region in M31. Only in the lower panel, 
with the more metal-poor clusters, do we see a signature 
like that of a kinematically hot classical bulge: there are 
roughly equal numbers of clusters in each quadrant, and 
we see that the velocity dispersion rises sharply close to 
the center, as we would expect for a centrally concen- 
trated classical bulge. However, as we shall show below, 
the starlight in this region is dominated by old metal- 
rich stars of near solar abundance, so these metal-poor 
clusters are not tracing the dominant component here. 

To summarize: we see strong evidence from the kine- 
matics of the metal-rich old clusters ([Fe/H]>-0.6) for 
both disk and bar kinematics. A number of the clus- 
ters within 2 kpc of the center of M31 show the kine- 
matic signature of xi orbits in a barred potential. The 
rest of these clusters (plus the other metal-rich clusters 
within 2 kpc of the major axis) show the cold kinematics 
of the disk. These k inema tics strongly confirm the re- 
sult o f [Beaton et all (|2007| ) and lAthanassoula k, Beatonl 
(|2006| ) that M31 has a bar whose inner parts constitute 
the boxy bulge which dominates its light in the inner few 
kpc. To our knowledge, this is the first clear detection 
of globular clusters with bar kinematics in any galaxy. 
However, there is one massive cluster (the Arches clus- 
ter) in the Milky Way which has a large space velocity 
(232 km s" 1 ) and is currently at a p rojected distance o f 
only 26 pc from the galactic center (jStolte et al.l 120081) . 
These authors note that the cluster could be on a transi- 
tional trajectory between x\ and X2 orbits in the Milky 
Way's barred potential, and may have been formed in a 
starburst triggered when a massive molecular cloud "col- 
lided on the boundary between x\ and xi orbits in the 
inner bar" . 

3.1. Relations between clusters and bulge/disk field stars 

We now consider the relationship between field stars 
and globula r clus ters in the inner regions of M31. 
iTrager et afl ([20001 ) studied the integrated light of M31's 
bulge, in a circular aperture of diameter 250 pc. They 
found a mean metallicity of +0 .2 dex, and a mean age of 
around 6 Gyr. More recently, iSaglia et al.l ([20101 ) have 
made a detailed study of the M31 bulge region using a 
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number of long-slit exposures with the HET. They find 
a mean metallicity around solar, and an age of around 
12 Gyr in the inner 1-2 kpc. (Note that they do see 
a metallicity gradient, reaching up to [Z/H]=+0.4, over 
the inner 200 pc, the region dominated by the classical 

b ulge.) 

ISaraiedini fc Jablonkal d200l used HST/WFPC2 ob- 
servations to produce a color-magnitude diagram for 
M31's bulge at 1.6 kpc from its center, and inferred 
a metallicity dis tribution which peaked near solar. 
lOlsen et al.l (|2006| ) summarized near IR color-magnitude 
diagrams from high spatial resolution studies of M31 to 
find that the stellar population in the inner few kpc was 
dominated by old, nearly solar-mctallicity stars. Inter- 
estingly, by comparing fields in the bulge with an inner 
disk field, they found no evidence for an age difference 
between bulge and disk. This is unsurprising if M31's 
bulge is dominated by a bar, since bar stars are merely 
inner disk stars which have become part of the bar pat- 
tern. 

The mean metallicity of the integrated light from field 
stars thus exceeds the mean metallicity of the globulars 
in the inner few kpc; it is closer to the mean of those 
with [Fe/H]>-0.6, which show either disk or bar kine- 
matics. (It has been suggested before that globular clus- 
ters are formed less efficiently in metal-rich populations: 
IStrader et al.l (|2005l ) calculate that the efficiencies differ 
by more than a factor of 10 in the Milky Way, by compar- 
ing metal-rich globular clusters to the bulge luminosity 
and metal-poor numbers to the halo luminosity. This 
number will not be changed radically if we substitute 
the thick disk luminosity for the bulge luminosity in this 
calculation.) 

Thus a simple picture can explain the existence of the 
metal-rich globular clusters in M31: they merely partic- 
ipated in the early formation of the inner disk and the 
onset of the bar instability. 
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We have discussed accurate kinematical data for old 
M31 clusters in its inner regions within 2 kpc of its ma- 
jor axis. The majority of the metal-rich clusters (those 
with [Fe/H] greater than -0.6) show disk kinematics, and 
many of the clusters within the innermost bar region have 
the signature of the X2 family. This clearly shows the ex- 
istence of an ILR and, to our knowledge, this is the first 
time it has been clearly shown usi ng stellar kinematics . 
In the only other known example, iTeuben et al.l (|1986f ) 
showed this using gas kinematics in the strongly barred 
galaxy NGC 1365. Our result also gives an estimate of 
the ILR location, which provides useful constraints for 
future dynamical studies of M31 since it could be used 
to set limits to the bar pattern speed. These metal-rich 
clusters share the population properties (metallicity and 
age) of the integrated light in the inner few kpc, which 
has been studied both via spectroscopy and via deep 
color-magnitude diagrams from HST and adaptive op- 
tics imaging. By contrast, clusters with [Fe/H] less than 
-0.6 within 2 kpc of the major axis show little rotational 
support and a velocity dispersion which increases as ra- 
dial distance to the center decreases. 

Our data do not probe the small region (200 pc) oc- 
cupied by M31's c lassical bulge in the description of 
iBeaton et all (|2007{) . so we cannot comment on its kine- 
matics. However, we caution against simply interpret- 
ing a high velocity dispersion in M31's inner few kpc 
as a bulge velocity dispersion and then using it to con- 
strain M31's black hole mass (as done most recently by 
iSaglia et al.l l2010h : the contribution of the bar, which 
dominates the light there, needs to be assessed. 
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